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UNDERSTANDING THE EQUILIBRIUM STRUCTURE OF 

CDM HALOS 
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Hugo Martel 3 



Abstract. N-body simulations find a universal structure for the ha- 
los which result from the nonlinear growth of Gaussian-random-noise 
density fluctuations in the CDM universe. This talk summarized our 
attempts to derive and explain this universal structure by analytical ap- 
proximation and simplified models. As an example, we show here that a 
fD spherical infall model involving a fluid approximation derived from 
the Boltzmann equation can explain not only the halo density profile 
but its phase-space density profile, as well. 



Q-f 1 Introduction 

o ■ , n 

Numerical N-body simulations of structure formation from Gaussian-random- noise 
CZ) , initial conditions in the CDM universe find a universal structure for halos. This 
' universality is a fundamental prediction of the CDM model, but our knowledge 
J> . is limited to the "empirical" N-body results, with little analytical understanding. 
In his talk, Shapiro summarized our attempts to fill this gap by a hierarchy of 
approximations, each simpler than the last: 1. 3D gas/N-body simulations of halo 
' formation from simplified initial conditions; 2. ID, spherical, analytical models 
using a fluid dynamics approximation derived from the Boltzmann equation; 3. 
An analytical hydrostatic equilibrium model which follows from the virialization 
of top- hat density perturbations. 

Most of the work described in that talk is s u mmar ized in Shapiro et al. (2004) 
and references therein and in lAhn fc Shapiro! l|2005|) , with the exception of our 
new results on the halo phase-space density profile, which we shall present here 
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for the first time. Due to length limitations, we shall limit this paper to just a 
few items from category 2 above. A more complete version of Shapiro's talk is 
available at the meeting website 1 . 

2 Universal Structure of CDM Halos: N-body Results 

CDM N-body halos show universal mass profiles. The same density profile fits 
halos from dwarf galaxies to clusters, independent of halo mass, of the shape 
of the density fluctuation power spectrum P(k), and of background cosmology: 
p(r)/p_2 = fcn(r/r- 2 ), where r_ 2 = radius where din p/d\nr = —2 and p_ 2 = 
p{ r -2) (e.g. Navarro et al. 2004) 2 . As r — > oo, p — > r~ 3 , while as r — > 0, p —> r~ a , 
1 < a < 1.5 (e.g. Navarro, Frenk, & White 1997; Moore et al. 1999). Diemand, 
Moore & Stadel (2004) report that 

PaSil ^ [r/r^[l + {r/r s rr-^ ^ 

with (a, P, 7) = (1, 3, 7) summarizes the fits to current simulations, with 7best-fit = 
1.16 ±0.14. 

The profiles of individual halos evolve with time. The halo mass grows as 
M(a) = Mooexpl-Saf/a], where 5" = [dlnM/dlna] (a = a f ) = 2 (Wechsler et 
al. 2002). The density profile concentration parameter, c = r 2 ooAs, also grows, 
according to c(a) = c(a,f)(a/af) for a > af (Bullock et al. 2001; Wechsler et al. 
2002), starting from c(o) < 3 — 4 for a < a/ (initial phase of most rapid mass 
assembly) (Tasitsiomi et al. 2004). 

CDM N-body halos show surprisingly isotropic velocity distributions. Halos 
have universal velocity anisotropy profiles, /3(r/r 2 oo), where (i(r) = l — (v^}/(2{v^}), 
with (3 = (isotropic) at r = 0, gradually rising to (3 ~ 0.3 at r = r 2 oo (e.g. Carl- 
berg et al. 1997). 

CDM N-body halos show universal phase-space density profiles. N-body results 
find p/<7y oc r _a P a , where a ps = 1.875 (Taylor & Navarro 2001), a ps = 1.95 (Rasia, 
Tormen and Moscardini 2004), and a ps = 1.9 ±0.05 (Ascasibar et al. 2004). (Also 
related: P(f) oc /-2-5±o.05. Arad; Dcke i ; & Klypin 2004). 

3 The Fluid Approximation: ID Halo Formation From Cosmological 
Infall 

The collisionless Boltzmann equation and the Poisson equation can be used to de- 
rive exact dynamical equations for CDM which are identical to the fluid conserva- 
tion equations for an ideal gas with adiabatic index 7 = 5/3, if we assume spherical 



1 http://www2.iap.fr/users/gam/yappa-ng/indcx.php?album=%2FIAP05%2F&image=Shapiro.pdf 
2 The weak mass-dependence suggested by Ricotti (2003) is an exception which remains to be 
confirmed. 
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Fig. 1. Self-similar Spherical Infall with e = 1/6. (a) (left) (top) Halo mass density 
versus radius for analytical, similarity solution in the fluid approximation, compared 
with best-fitting NFW, Moore and a/37 profiles with (a,/3,7)=(l,3,7). (bottom) fractional 
deviation of fits from self-similar solution pss- (b) (right) Halo phase-space density versus 
radius for analytical similarity solution compared with best fitting power-law r -1 ' 91 . 



symmetry and a velocity distribution which is both skewless and isotropic, assump- 
tions which approximate the N-body results reasonably well (Ahn & Shapiro 2005). 
We have used this fluid approximation to show that most of the universal proper- 
ties of CDM N-body halos described above can be understood as the dynamical 
outcome of continuous cosmological infall. 

3.1 Self-similar gravitational collapse: the spherical infall model 

In an EdS universe, scale- free, spherically symmetric perturbations 5M/M oc M~ c 
result in self-similar structure formation. Each spherical mass shell around the 
center expands until it reaches a maximum radius r ta and recollapses, r ta cx t^, 
£ = (6e + 2)/(9e). There are no characteristic length or time scales besides r ta 
and Hubble time t. For cold, unperturbed matter, this results in highly supersonic 
infall, terminated by a strong, accretion shock which thermalizes the kinetic energy 
of collapse: r s h oc k(i) °c ^ta(i)- The spherical region bounded by the shock is 
roughly in hydrostatic equilibrium, a good model for virialized halos. 

Consider halo formation around peaks of the Gaussian random noise primordial 
density fluctuations. If P{k) oc k n , then ^cr-peaks with v > 3 have simple power- 
law profiles for accumulated overdensity inside r, A (r) = SM/M oc r~("+ 3 ) oc 
^j-(n+3)/3 ( e _g Hoffman & Shaham 1985), which implies self-similar infall with 
e = (n + 3)/3. For ACDM, galactic halos are well approximated by n = —2.5 ±0.2 
for 10 3 < M/Mq < 10 11 . By applying the fluid approximation to the problem 
of self-similar spherical infall with e = 1/6 (i.e. n ff = —2.5), Ahn & Shapiro 
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Fig. 2. Non-Self-Similar Spherical Infall. (a) (left) (top) Halo mass density versus radius 
at epoch a/at = 3.85, according to fluid approximation solution for the non-self-similar 
spherical infall rate which makes halo mass grow in time like Wechsler et al. (2002) fitting 
formula, compared with best-fitting NFW and Moore profiles; (bottom) Corresponding 
halo rotation curves, (b) (right) Concentration parameter versus scale factor for the best- 
fitting NFW profiles at each time-slice during the evolution of the fluid approximation 
solution for time-varying spherical infall at the Wechsler et al. (2002) rate. 

(2005) derived a ID, analytical solution for halo formation and evolution, in which 
''shock cx r ta oc i 2 , and M a ( 4 . The resulting self-similar halo density profile 
inside the accretion shock agrees with that of CDM N-body halos, with a best- 
fit a/?7-profile which has (a, (3, 7) = (1,3,1.3) (see Figure 1(a)). As we show in 
Figure 1(b), this analytical similarity solution for e = 1/6 also derives the universal 
phase-space density profile found for CDM N-body halos, p/o\ oc r~ aps , a ps ~ 1.9. 

3.2 Non-self-similar infall: Mass Assembly History and the Evolution of CDM 
N-body Halo Profiles 

Self-similar infall may provide a good explanation for some halo properties, but it 
cannot explain how profile shapes change with time and depart from self-similarity. 
To do this, we have derived the perturbation profile that leads to the non-sclf- 
similar halo mass growth rate reported for CDM N-body halos by Wechsler et al. 
(2002) and used the fluid approximation to derive the halo properties that result 
(Alvarez, Ahn, & Shapiro 2003). We solved the fluid approximation equations 
by a ID, spherical, Lagrangian hydro code. These solutions explain most of the 
empirical CDM N-body halo properties and their evolution described above as 
a dynamical consequence of this time-varying, but smooth and continuous infall 
rate. The halo density profiles which result are well fit by (and intermediate be- 
tween) NFW and Moore profiles (over the range r/r 2 oo > 0.01) (see Figure 2(a)). 
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Fig. 3. Phase-space density profiles for the same non-self-similar infall solution plotted 
in Figure 2, for a/at = 1, 2.1, and 6.1, with arrows showing locations of V2oo at each 
epoch, along with best-fitting power-law r -1 ' 93 . 



These halo density profiles evolve just like CDM N-body halos, too. The halo 
concentration parameter grows with time just like CDM N-body halos (see Figure 
2(b)). In addition, these solutions yield a halo phase-space density profile, p/cr^, 
in remarkable agreement at all times with the universal profile reported for CDM 
N-body halos (see Figure 3). We therefore conclude that the time- varying mass 
accretion rate, or equivalently the shape of the initial density perturbation, is the 
dominant influence on the structure of CDM halos, which can be understood sim- 
ply in the context of spherical collapse and the accretion of a smoothly-distributed, 
spherically-symmetric collisionless fluid. 

This work was supported NASA Astrophysical Theory Program grants NAG5- 
10825, NAG5-10826, NNG04GI77G, Texas Advanced Research Program grant 
3658-0624-1999, and a Department of Energy Computational Science Graduate 
Fellowship to M.A.A. 
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